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INTRODUCTION 


I 

I 


: 

: 


The  efforts  described  herein  are  part  of  a project  to  develop  and  im- 
plement a data  processing  software  system  for  the  management  and 
analysis  of  data  received  from  the  AFGL  MESA  (Miniature  Electro- 
static Accelerometer)  experiment  flown  aboard  the  NASA  Atmosphere 
Explorer  (AE)  series  of  satellites.  To  accomplish  this  task  a Xerox 
Sigma-9  computer  system  — dedicated  to  the  AE  project  — is  being 
utilized  in  a time-shared  environment.  Telemetry  data  from  the 
major  experiments  aboard  AE  are  transmitted  from  remote  stations 
to  NASA’s  Goddard  Space  Flight  Center  (GSFC),  and  these  data  are 
then  stored  on  mass  storage  devices  for  use  by  experimenter's  soft- 
ware. Reduced  and  analyzed  experiment  data  are  later  stored  on 
mass  storage  devices  for  use  by  all  other  experimenters  and  theo- 
retical analysts. 

The  MESA  Data  Reduction  System  (DRS)  which  has  been  developed  is 
capable  of  extracting  raw  telemetry  data  from  the  AE  data  base,  edit- 
ing and  temperature-correcting  the  telemetry  data,  extracting 
atmospheric  drag  values  utilizing  digital  filtering  techniques,  and 
calculating  atmospheric  density  and  wind  data.  In  addition,  the  MESA 
DRS  calculates  Jacchia  71  model  density  values,  lists  and  displays 
calculated  parameters  on  printer  plots  and/ or  microfilm,  and  stores 
output  data  in  files  for  on-line  use  by  other  experimenters.  This  sys- 
tem is  presently  being  utilized  in  a production  environment  to  produce 
reduced  MESA  density  data  from  three  AE  satellites  on  a timely  basis. 

This  report  will  describe  the  AE  satellite,  the  AE  data  system,  the 
MESA  experiment,  and  the  MESA  data  processing  system. 
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1. 1 The  AE  Satellites 

The  Atmosphere  Explorer  satellites  are  a series  of  maneuverable, 
unmanned  spacecrafts  whose  purpose  is  to  investigate  the  physical 
properties,  dynamics  and  photochemical  processes  in  the  upper 
atmosphere.  The  spacecrafts  differ  from  the  usual  unmanned 
scientific  satellites  in  that  each  contains  an  on-board  propulsion 
system  which  permits  variation  of  perigee  and  apogee  altitudes;  in 
the  team  approach  taken  by  investigators  to  analyze  and  compare 
data;  in  the  normal  AE  spacecraft  and  data-collecting  operations; 
and  in  the  concommitant  rapidity  with  which  data  must  be  forwarded, 
processed  and  analyzed. 


Linked  through  a dedicated  ground-based  computer  operated  in  a 
timeshared  environment  with  scientists  in  various  parts  of  the  coun- 
try, the  AE  satellites  carry  as  many  as  fourteen  scientific  experi- 
ments designed  to  provide  data  on  natural  processes  and  phenomena 
in  the  earth's  thermosphere. 

The  AE-C  spacecraft  was  launched  from  Vandenburg  Air  Force  Base, 
California  on  16  December  1973  at  06:18  GMT  into  an  elliptical  orbit 
with  apogee  at  about  4000  km,  perigee  at  156  km,  and  an  inclination 
of  68.4  degrees.  AE-D  was  launched  from  Vandenberg  Air  Force 
Base  on  6 October  1975  at  09:00:49  GMT  into  a polar  orbit.  Initial 
apogee  was  about  3800  km,  perigee  156. 4 km  and  the  inclination  was 
90. 1 degrees.  The  third  satellite  in  the  series,  AE-E,  was  launched 
from  Cape  Canaveral,  Florida  on  20  November  1975  at  02:06:48  GMT. 

It  was  injected  into  an  elliptical  equitorial  orbit  with  apogee  at  3025  km, 
perigee  156. 9 km,  and  an  inclination  of  19. 7 degrees. 

The  AE  mission  objective,  to  study  phenomena  in  the  atmosphere  at 
altitudes  above  120  km,  was  to  be  accomplished  by  each  satellite  in 
two  orbital  phases:  elliptical  orbit  and  circular  orbit.  During  the 


first  few  months  in  orbit,  the  perigee  of  each  satellite  was  to  be  peri- 
odically lowered  to  within  the  altitude  ranges  of  130  km  to  160  km. 
This  was  to  provide  data  on  the  thermosphere  — a region  that  has  been 
essentially  unexplored  on  a continuous  basis.  After  these  excursions 
the  spacecrafts  were  to  be  moved  slowly  into  circular  orbits  at  speci- 
fied altitudes. 

AE-C  successfully  completed  its  elliptical  phase  and  was  circularized 
in  December  1975,  at  about  240  km.  AE-D  and  AE-E  are  presently  in 
the  elliptical  phase  of  their  respective  missions. 

This  combination  of  both  elliptical  and  circular  orbits  taken  on  polar 
and  equitorial  orbiting  satellites  will  provide  global  coverage  of  the 
upper  atmosphere.  Hopefully,  this  will  aid  in  separating  such  atmos- 
pheric effects  due  to  altitude,  latitude,  and  solar  radiation. 

Reduced  geophysical  data  are  stored  in  the  central  computer  and  are 
shared  among  all  AE  scientists.  This  permits  unified  attacks  on  the 
relevant  problems  and  also  allows  a minimal  delay  in  the  publication 
of  results. 

A more  detailed  description  of  the  AE  satellite  may  be  found  in  Ref- 
erence [1J. 
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1 . 2 The  AE  Data  System 


The  AE  Central  Computer  located  at  Goddard  Space  Flight  Center 
(GSFC),  Greenbelt,  Maryland,  consists  of  a Xerox  Sigma-9  computer 
system  under  the  Control  Program  No.  5 (CP-5)  Operating  System. 
The  Sigma-9  is  equipped  with  1. 28  megabytes  (320K  32-bit  words)  of 
core  storage,  10. 6 megabytes  of  Rapid  Access  Device  (RAD)  on-line 
storage,  800  megabytes  of  high  speed  disk  on-line  storage,  and  10 
9-track  drives  (seven  1600-bpi  drives  and  three  800-bpi  drives).  A 
Xerox  Sigma-3  computer,  which  is  connected  to  the  Sigma-9  com- 
puter, is  used  as  a controller  to  interface  with  various  Sanders  810 
programmable  terminal  systems  located  at  GSFC  and  the  experiment- 
ers remote  sites.  Hardwire  data  links  to  the  Sigma-9  provide  space- 
craft telemetry  data  input  from  the  Input  Processor  Computer  (Xerox 
Sigma-5)  and  orbit/attitude  data  input  from  the  Orbit/ Attitude  Proces- 
sor Computer  (EBM  S/360-95). 

The  AE  Central  Computer  is  dedicated  to  the  support  of  the  AE-C, 
AE-D,  and  AE-E  missions.  In  particular,  most  of  the  processing 
on  the  Central  Computer  involves  the  reduction  and  analysis  of  data 
from  on-board  experimental  equipment.  The  development  of  software 
systems  to  perform  these  functions  is  the  responsibility  of  individual 
investigators.  To  simplify  this  development  and  to  ensure  flexibility 
in  data  handling  and  sharing  between  investigators,  standardized  data 
bases  have  been  established  for  the  Central  Computer  which  are  ac- 
cessible through  the  Fortran  language  by  using  custom  designed 
interface  routines.  The  data  bases  utilized  by  Investigators  are: 

a)  Time-smoothed  telemetry  data 

b)  Orbit/attitude  data 

c)  Magnetic  solar  activity  data 

d)  Geophysical  unit  data  (separate  file  structures  for  each 
experiment  after  initial  data  reduction) 


e)  Unified  abstract  data  (summary  results  from  each  ex- 
periment combined  in  a single  file  structure). 

These  data  bases  are  maintained  on  disk  and  tape  storage.  In  general, 
disk  storage  is  used  for  the  most  recent  data  available,  and  tape  stor- 
age is  used  for  permanent  storage  as  well  as  backup  for  the  data  on 
disk.  Figure  1 illustrates  the  location  of  disk  storage  to  those  data 
bases  requiring  large  amounts  of  space. 

A variety  of  software  resources  are  available  on  the  Central  Computer 
which  simplify  monitoring  as  well  as  processing  these  data  bases. 

The  majority  of  these  resources  are  included  in  the  Data  Management 
Facility  (DMF)  developed  under  the  auspices  of  the  GSFC  Information 
Processing  Division.  These  responsibilities  include: 

a)  Original  creation  and  maintenance  (updating)  of  the  data 
bases. 

b)  Management  of  data  base  storage  space  (disk  and  tape). 

c)  Software  interfaces  (FORTRAN  callable)  for  reading  and 
(when  applicable)  writing  the  data  bases. 

d)  On-line  query  capabilities  for  data  base  status  and  sum- 
mary data  information. 

e)  Utility  programs  for  data  base  manipulation. 

A more  detailed  description  of  the  Data  Management  Facility  and  its 
use  may  be  found  in  Reference  [2], 
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DISK  STORAGE 


SYSTEM  USE 

150  MEGABYTES  ALLOCATION 

^ALL  TELEMETRY  DATA  FROM  THE  SIX  MOST  RECENT  ORBITS., 
50  MEGABYTES  ALLOCATION 

ORBIT  ATTITUDE  INFORMATION  FOR  UP  TO  1 YEAR 

50  MEGABYTES  ALLOCATION 

UNIFIED  ABSTRACT  DATA  FILE  FOR  UP  TO  5 MONTHS  OF 
HEADER  AND  OTHER  DATA  FROM  ALL  INVESTIGATORS 
ON  15  SECONDS  TIME  INCREMENTS.  ASSUM  ING  30  PER 
CENT  DUTY  CYCLE 

100  MEGABYTES  ALLOCATION 


16  MEGABYTES  OF  PERMANENT  STORAGE  ASSIGNED  TO 
EACH  OF  1 7 INVESTIGATORS  PLUS  45  MEGABYTES  TO  BE 
TEMPORARILY  ASSIGNED  AMONG  INVESTIGATORS  FOR 
SPECIFIC  PROBLEMS 

_300  MEGABYTES  ALLOCATION 

DURING  PRIME  COMPUTER  SHIFT.  AREA  WILL  BE  OPEN 
SO  THAT  INDIVIDUAL  INVESTIGATORS  CAN  CHOOSE  ANY 
TELEMETRY  DATA  TO  BE  LOADED  FOR  ACCESS.  DURING 
EVENING  SHIFTS  AREA  WILL  BE  ASSIGNED  TO  SPECIFIC 
24-HOUR  8L  OCXS  OF  TELEMETRY  TO  BE  PROCESSED  BY 
INVESTIGATOR  PRODUCTION  BATCH  PROGRAMS.  TWO 
SUCH  BLOCKS  CAN  BE  PROCESSE  D CONCURR  ENTLY 

150  MEGABYTES  ALLOCATION 


Figure  1.  AE  Central  Computer  Disk  Storage  Allocation 
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1.3  The  MESA  Experiment 

The  MESA  experiments  on  the  AE  satellites  were  designed  to  provide 
accurate  neutral  atmospheric  density  measurements  by  measuring 
satellite  deceleration  caused  by  aerodynamic  drag.  The  individual 
MESA  unit  is  a uniaxial  accelerometer.  It  consists  of  an  electro- 
statitically  suspended  proof  mass  which  is  also  electrostatically 
rebalanced  along  a sensitive  axis.  The  MESA  determines  an  applied 
acceleration  from  the  electrostatic  force  required  to  recenter  the 
proof  mass. 

The  output  of  the  MESA  is  a digital  pulse  rate  proportional  to  the 
sensed  input  acceleration.  Vehicle  dynamics,  the  momentum  wheel, 
propulsion  system  thrusting,  and  instrument  motions  provide  input 
"noise"  accelerations.  These  noise  accelerations  are  to  be  removed 
in  the  data  analysis  in  order  to  retrieve  the  desired  "signal"  accelera- 
tions due  to  atmospheric  drag. 

Three  single-axis  miniature  electrostatic  accelerometers  were  flown 
on  each  of  the  AE  missions.  Two  sensors  were  placed  with  their 
sensitive  axes  in  the  spacecraft  xy  plane.  Sensor  yx  was  placed  at 
the  center  of  gravity  with  the  sensitive  axis  45°  from  the  +x-axis  in 
the  +x,  -y  quadrant,  and  sensor  xy  is  placed  adjacent  to  yx  with  the 
sensitive  axis  45°  from  the  +x-axis  in  the  +x,  +y  quadrant  and  with  the 
center  of  the  proof  mass  on  the  spin  axis.  These  two  instruments  are 
used  to  provide  redundant  capability  for  density  measurement,  partic- 
ularly during  the  despun  mode,  as  well  as  monitoring  the  orbit-adjust 
propulsion  system.  The  third  sensor  was  placed  with  its  sensitive 
axis  along  the  spacecraft  z-axis.  This  sensor  can  be  used  to  monitor 
spacecraft  roll,  particularly  in  the  despun  mode  of  operation. 


! 


The  MESA  sensors  have  the  capability  of  being  commanded  into  any 
of  three  sensitivity  ranges: 


■ 
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Range 

A 

B 

C 


Full  scale  (g's) 

8 x l(f  3 
4 x 10~4 
2 x l(f 5 


Range  A was  commanded  for  satellite  thrust  monitoring  and  was  not 
used  for  density  measurement.  Emphasis  in  this  report  will  be  given 
to  data  taken  from  sensors  in  range  B,  since  most  density  data  were 
taken  in  this  range.  The  altitude  region  in  which  drag  can  be  meas- 
ured for  each  sensitivity  range  is  as  follows: 


Range  A 
Range  B 
Range  C 


120  to  160  km 
122  to  280  km 
180  to  320  km 


Each  sensor  determines  the  component  of  drag  along  a single  axis 
only. 


Once  the  component  of  the  drag  acceleration  is  calculated  for  each 
single  axis,  the  total  drag  acceleration  may  be  calculated  in  one  of 
three  ways.  First,  if  three  instrument  axes  are  active,  the  total 
drag  acceleration  A^  is 


where 


ai 

8, 


drag  acceleration  component  of  the  i^1  accelerometer 
i^  accelerometer  scale  factor. 


Secondly,  if  any  one  accelerometer  is  turned  on,  the  total  drag  acceler- 
ation may  also  be  derived  from  that  single  axis  separately  by 
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A = * 

D s.  cos  Of^  * 

where  cos  a.  is  the  angle  between  the  accelerometer  sensitive 
axis  and  the  satellite's  velocity  vector. 

It  is  important  to  note  that  since  the  accelerometer  outputs  are  aver- 
aged accelerations  sensed  over  the  sample  interval  (one-quarter  sec- 
ond), the  cos  a.  here  should  be  an  averaged  value  taken  over  the 
sample  interval. 


In  the  third  case,  if  the  xy  and  yx  accelerometers  are  active,  the  total 
drag  acceleration  may  be  calculated  by  locating  those  times  when  the 
z-accelerometer’s  sensitive  axis  is  making  a 90°  attack  angle  with 
the  satellite's  velocity  vector.  In  that  case,  the  drag  along  that  axis 
should  be  zero  and  the  total  drag  acceleration  can  be  calculated  by 


*■>{  z (S 

\i=xy,yx 


In  any  case,  density  is  then  derived  from  the  following: 


CDaV 

where  p atmospheric  density 

m/a  ■ satellite  mass  to  area  ratio 
Ad  - total  drag  acceleration 
CQ  - satellite  drag  coefficient 
V - satellite  velocity. 


A complete  description  of  the  MESA  experiment  may  be  found  in 
Reference  [3], 
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1.3.1  Digital  filtering  analysis 

The  approach  taken  to  extract  the  message  signal  accelerations  from 
the  total  signal  accelerations  is  provided  by  statistical  communication 
theory. 

1. 3. 1. 1 Despun  orbits  - In  the  case  when  the  satellite  is  despun,  at- 
mospheric drag  information  assuming  no  wave  motions,  are  sensed 
as  zero  frequency  (DC)  accelerations  whose  amplitudes  vary  between 

Q A 

10  and  5 x 10  g's  depending  upon  satellite  altitude.  Satellite 

nutation  modifies  the  MESA  output  signal  with  a sinusoidal  modulation 

whose  frequency  is  0. 14-0.2  Hz  and  whose  amplitude  is  about 
-4 

1.8  x 10  g's  per  degree  of  nutation  of  the  XY  sensor.  Typically, 
nutation  accelerations  vary  between  10  and  10  g's. 

Sensor  outputs  are  further  modified  by  noise  due  to  the  rotation  of  the 

satellite's  momentum  wheel  assembly  (MWA).  The  momentum  wheel 

rotation  causes  noise  accelerations  at  high  frequencies  which  are 

sensed  by  MESA.  However,  due  to  the  telemetry  data  sampling  rate, 

the  modulations  for  the  most  part  appear  at  1.  5—2.0  Hz.  This  effect 

is  known  as  aliasing  and  is  discussed  in  Reference  [2].  The  amplitude 

-6  . 

of  the  MWA  noise  accelerations  is  about  8 x 10  g's. 

Instrument  bias  causes  a DC-output  component  whose  amplitude  is 
—6 

about  1.4  x 10  g's.  In  the  despun  case,  centripetal  accelerations 
sensed  by  MESA  are  negligible. 

Figure  2 displays  a small  segment  of  raw  MESA  data  from  the  XY 
sensor  taken  on  orbit  2217  at  an  altitude  of  about  170  km.  Modula- 
tions on  the  output  signal  due  to  drag,  nutation,  MWA  noise  and 
instrument  bias  are  shown. 

Results  of  a power  spectrum  analysis  of  the  XY  sensor  data  are  dis- 
played in  Figure  3 showing  increased  power  due  to  (a)  atmospheric 


Figure  2.  Raw  MESA  Data  from  XY-Seneor 
Orbit  2217  Day  74180 
Despun  170  KM  Downleg 
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drag  at  DC,  and  (b)  momentum  wheel  and  nutation  noise  at  1. 95  Hz 
and  0. 156  Hz,  respectively.  Base  line  data  are  displayed  at  the  rela- 
tive white  noise  level. 

It  is  essential  that  the  filter  used  to  reduce  these  data  have  near  100% 
response  at  low  frequencies  and  minimum  response  at  those  frequen- 
cies where  nutation  and  momentum  wheel  noise  appear;  that  is,  that 
filter  side-lobes  be  kept  to  a minimum.  Figure  4 illustrates  the  filter 
chosen  to  accomplish  these  requirements.  The  filter  is  designated 
"15-10"  because  it  is  a low  pass  filter  whose  response  is  near  100% 
from  DC  to  15  seconds  and  is  minimal  from  10  seconds  to  .5  seconds. 

Since  the  "15  - 10"  filter  chosen  allows  DC  components  to  remain  in  the 
data,  what  remains  after  filtering  are:  (a)  drag  accelerations  and  (b) 
instrument  bias.  Bias  values  are  removed  by  considering  the  filtered 
output  only  in  those  regions  where  atmospheric -drag  is  negligible,  that 
is,  at  high  altitudes.  Once  bias  values  are  determined,  they  are  sub- 
tracted from  the  filtered  output  and  the  remaining  values  are  due  to 
atmospheric  drag.  Figure  5 displays  MESA  density  data  from  orbit 
2217  after  filtering  has  been  done. 

1. 3. 1. 2 Spinning  orbit  - We  now  consider  the  case  when  the  satellite 
is  spinning,  nominally  at  4 rpm.  In  this  configuration,  although  atmos- 
pheric drag  information  would  normally  appear  at  low  frequencies  (near 
DC)  the  satellite  spinning  causes  the  drag  measurements  to  be  "chopped" 
at  the  spinning  rate.  Mathematically,  this  is  equivalent  to  multiplying 
the  measured  drag  signal  by  a sine  wave  at  the  spin  frequency.  Equiv- 
alently, the  power  spectrum  of  the  drag  signal  is  shifted  in  frequency  to 
the  spin  frequency. 

-8  -4 

The  amplitude  of  the  drag  information  varies  between  10  and  5 x 10 
g's  depending  upon  the  satellite  altitude.  In  the  spinning  mode  the  spin 
axis  nutation  modifies  the  MESA  output  as  a 0. 1 - 0. 11  Hz  sinusoid.  As 
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Figure  3.  Power  Spectrum  Analysis 
Orbit  2217  Day  74180 
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Figure  4.  Filter  '15—10'  Response  Curve 
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in  the  despun  mode,  nutation  accelerations  average  about  1.8  x 10  g's, 
per  degree  of  nutation  on  the  MESA  XY  sensor. 

Momentum  wheel  related  accelerations  for  the  most  part  again  appear 

as  high  frequency  modulations  at  0. 6 - 0. 8 Hz  with  amplitudes  averaging 
— fi 

about  8 x 10  g's.  In  addition,  accelerations  due  to  movement  of  other 
sensors  sometimes  appear  at  frequencies  lower  than  the  spin  frequency, 
typically  about  0. 033  Hz. 

Figure  6 displays  typical  output  data  from  MESA  on  orbit  2437  when  the 
satellite  was  spinning  at  0.067  Hz.  Power  spectrum  analysis  results 
from  the  XY  sensor  data  are  shown  in  Figure  7.  Drag  information  is 
displayed  as  increased  power  at  the  spin  frequency,  with  nutation  effects 
and  momentum  wheel  noise  indicated. 

To  extract  drag  information  from  the  total  MESA  signal  in  the  spinning 
mode,  a "band-pass"  filter  was  designed  with  the  characteristics  shown 
in  Figure  8.  This  filter  removes  (a)  unwanted  nutation  accelerations 
which  appear  at  about  0. 102  Hz,  (b)  some  motions  of  other  instruments 
at  about  0. 033  Hz,  (c)  bias  and  centripetal  accelerations  appearing  near 
DC,  and  (d)  momentum  wheel  noise  at  0. 6 - 0. 8 Hz.  At  the  same  time 
it  is  centered  at  0. 067  Hz  to  allow  atmospheric  drag  information  to  pass. 

The  filter  parameters  were  chosen  to  ensure  that  the  filter  would  de- 
scribe the  maximum  variation  in  atmospheric  drag  by  passing  the  drag 
signal  information,  centered  at  the  satellite  spin  frequency,  having  a 
bandwidth  within  the  filter  bandwidth.  That  is,  the  drag  signal  bandwidth 
is  within  . 013  Hz  of  DC  before  modulation  by  the  spin  frequency. 

The  spectrum  of  this  drag  signal  after  spin  modulation  will  reside  within 
the  filter  bandwidth.  At  the  same  time  all  non-atmospheric  noise  accel- 
erations at  frequencies  of  less  than  0. 05  Hz  and  greater  than  0. 1 Hz  will 
be  removed  by  the  filter. 
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Figure  7.  Power  Spectrum  Analysis 
Orbit  2437  Day  74198 
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Since  DC  components  are  removed  from  the  data  by  the  bandpass  filter 
(as  indicated  by  the  response  curve  in  Figure  8)  instrument  bias  and 
centripetal  accelerations  are  filtered  out  and  hence  need  not  be  sepa- 
rately removed. 

Figure  9 displays  MESA  density  data  from  orbit  2437  after  filtering  has 
been  done.  A complete  description  of  the  filtering  techniques  utilized 
may  be  found  in  Reference  [4], 


FIGURE  9.  ATMOSPHERIC  DENSITY  DATA  ORBIT  2437  DAY  74198  SPINNING 
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2.  DATA  PROCESSING 

2. 1 MESA  Data  Reduction  System 

The  MESA  data  reduction  system  (DRS)  was  initially  written  (prior  to 
the  launch  of  AE-C)  as  an  overlay  system  of  computer  programs. 
Although  the  structure  has  not  changed,  the  programs  have  been  pe- 
riodically modified  to  their  present  forms  and  are  now  being  utilized 
in  a production  environment  to  produce  MESA  density  data  from  AE-C, 
AE-D,  and  AE-E  on  a timely  basis. 

Presently,  the  MESA  DRS  is  made  up  of  a file  naming  routine,  ROTATE, 
a setup/control  program,  DRIVER9,  and  six  primary  processing  rou- 
tines, which  perform  the  necessary  MESA  data  reduction  and  display 
functions.  The  primary  processing  routines  are:  RAWDATA,  DRAG, 
MODPAP,  FILM,  UAMAKE,  and  GUMAKE.  Descriptions  of  the  rou- 
tines which  comprise  the  MESA  DRS  are  given. 

2. 1. 1 ROTATE  — ROTATE  is  a small  preprocessing  routine  whose 
sole  purpose  is  to  generate  suffixes  which  will  be  appended  to  all  output 
file  names  of  a particular  DRS  run.  These  suffixes  guarantee  that  each 
file  will  have  a unique  name.  ROTATE  allows  for  either  (a)  the  same 
job  to  be  run  a series  of  times  without  any  loss  of  valuable  output,  or 
(b)  simultaneously  running  two  executions  of  the  same  job  with  different 
input  data. 

2.1.2  DRIVER9  - DRIVER9  is  a setup  and  control  program  which  gen- 
erates namelist  input  values  for  the  primary  MESA  processing  routines, 
utilizing  as  input  user-created  directive  records  called  Orbit  Request 
Cards  (ORCs).  Various  parameters  allow  full  control  over  each  primary 
program  in  the  DRS,  with  additional  control  over  secondary  processors 
to  promote,  demote,  and/or  remove  the  associated  GU-files  (Geophysical 
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Unit  files).  For  every  set  of  namelists  generated,  a corresponding  re- 
port is  written,  listing  any  errors  detected  and  any  actions  taken. 

There  are  two  types  of  ORC  parameters,  Telemetry  (TM)  and  Geo- 
physical Unit  (GU).  Telemetry  parameter  sets  are  used  to  process 
data  from  the  telemetry  data  file  and  usually  indicate  default  process- 
ing from  raw  telemetry  data  to  final  density  values.  The  form  of  the 
TM  parameter  syntax  is: 

SDATE,  STIME,  EDATE,  ETIME  [,  SATID) 

and  is  broken  down  as  follows: 

SDATE  5-digit  Julian  start  date  of  data  to  be  processed, 
form  YYDDD. 

STIME  8-digit  start  time  of  data,  in  milliseconds, 
range  0-86399999. 

EDATE  5-digit  Julian  stop  date  of  data.  If  00000,  assumed 
to  be  the  same  as  SDATE. 

ETIME  8-digit  stop  time  of  data,  in  milliseconds. 

SATID  1-character  satellite  identifier,  either  C,  D,  or  E. 

This  parameter  is  optional.  If  omitted,  the  default 
SATID  is  E. 

Geophysical  Unit  (GU)  parameter  sets  are  used  to  process  data  from 
MESA  GU-files,  that  is,  data  which  has  been  previously  processed  in 
some  form.  The  form  of  the  GU  parameter  syntax  is:  *GU  ("file 
name")  followed  by  any  options  and/or  special  requests  which  need  be 
made.  Default  processing  is  from  pre-processed  temperature  cor- 
rected MESA  data  to  final  density  values. 


Figure  10  illustrates  both  TM  and  GU  parameter  requests.  The  first 
two  requests  are  for  TM  data  to  be  processed  for  day  74198  (that  is, 
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17  July  1974)  with  the  GUMAKE,  UAMAKE,  and  FILM  primary  routines 
to  be  executed  as  well  as  the  default  routines  RAWDATA,  DRAG,  and 
MODPAP. 


Requests  3 and  4 are  for  GU  data  to  be  processed  from  the  files  speci- 
fied. For  all  four  requests  namelist  input  values  will  be  generated  for 
the  primary  processing  routines  to  read,  and  descriptions  of  the  orbits 
which  these  requests  will  process  are  given. 

The  remaining  routines  are  the  primary  processing  routines  of  the  MESA 
DRS.  Following  is  a list  of  each  routine  and  the  types  of  input  and  output 
files  - telemetry  (TM)  geophysical  unit  (GU)  and  unified  abstract  (UA)  - 
which  each  requires: 


Name 

Input 

Output 

RAWDATA 

TM 

GU 

DRAG 

GU 

GU 

MODPAP 

GU 

GU 

FILM 

GU 

film  plots 

UAMAKE 

GU 

UA 

GUMAKE 

GU 

GU 

Following  is  a brief  description  of  each  routine. 

2.1.3  RAWDATA  - RAWDATA  reads  the  telemetry  file,  extracts  all 
MESA  information  for  those  times  when  the  satellite  is  below  700  kilo- 
meters, determines  on/off  and  sensitivity  status,  constructs  MESA 
sensor  outputs  from  the  TM  words,  determines  accelerometer  temper- 
atures and  temperature-corrects  measured  data,  edits  erroneous  data 
values,  replaces  missing  data  values,  calculates  appropriate  orbit/ 
attitude  information,  and  outputs  these  data  to  a GU  file  for  use  by  other 
primary  programs. 


(' 

I 
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Printer  output  from  RAWDATA  includes  accelerometer  temperature  in- 
formation, satellite  and  momentum  wheel  spin  rates,  and  various 
statistical  values  describing  the  quality  of  the  TM  data  processed  by 
this  program. 

Figure  11  is  an  example  of  the  two  page  printer  report  generated  by 
RAWDATA. 

2. 1.4  DRAG  — This  routine  is  the  major  processing  program  of  the 
DRS,  and  it  requires  as  input  the  GU  file  created  by  RAWDATA.  DRAG 
extracts  atmospheric  drag  information  from  the  temperature-corrected 
data  by  utilizing  various  digital  filters  and/or  data  modelling  techniques, 
calculates  atmospheric  density  values  by  one  of  the  methods  described 
in  Section  1.3,  determines  any  transverse  wind  components,  and  stores 
this  information  back  on  the  GU  file. 

In  addition,  DRAG  has  the  capability  of  calculating  attitude  information 
utilizing  the  MESA  output  data.  This  is  done  periodically  to  test  the 
accuracy  of  the  spacecraft  attitude  system. 

Printer  output  from  DRAG  includes  calculated  instrument  bias  values, 
MESA  on/off  and  range  status  and  coefficients  of  a least  squares  fit 
done  to  the  calculated  density  data.  These  coefficients  will  later  be 
used  to  generate  a printer  plot  display  of  the  density  versus  altitude 
profile  for  the  particular  orbit  processed. 

< Figure  12  illustrates  the  one  page  printer  report  generated  by  DRAG. 


2. 1.  5 MODPAP  — MODPAP  requires  the  GU  file  created  by  DRAG  as 
input  to  perform  its  three  main  functions.  It: 

(a)  calculates  and  stores  (in  the  GU  file)  Jacchia  71  model 
atmpsphere  density  values, 


tO*T*  Misr 
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(b)  provides  various  printer  plot  displays  of  MESA  generated 
atmospheric  parameters, 

(c)  generates  a comprehensive  listing  of  particular  orbit,  atti- 
tude, and  MESA  generated  parameters. 

Figure  13  displays  the  report  generated  by  the  MODEL  section  of 
MODPAP.  An  example  of  the  printer  plot  density  profile  which  may  be 
generated  by  MODPAP  is  given  in  Figure  14.  Figure  2 (Section  1. 3) 
is  an  example  of  a raw  data  printer  plot  displayed  as  a function  of  GMT. 
Figure  15  (two  pages)  illustrates  the  type  of  report  generated  by  the 
GU-list  section  of  MODPAP. 

2.1.6  UAMAKE  - UAMAKE  reads  the  density  values  generated  by 
DRAG  from  the  GU-file  and  writes  density  and/or  wind  values  into  their 
proper  position  in  the  Unified  Abstract  (UA)  file  for  use  by  other  experi- 
menters. Since  UA  data  is  stored  at  exact  15-second  intervals  by  all 
experiments,  UAMAKE  performs  interpolating  functions  to  generated 
density  values  at  the  required  times. 

Printer  output  from  UAMAKE  is  a listing  of  MESA  generated  density 
values  at  the  required  UA  times. 

Figure  16  displays  the  printer  report  generated  by  the  UAMAKE  routine. 

2.1.7  GUMAKE  - The  main  function  of  GUMAKE  is  to  compact  the  GU 
file  created  by  the  other  primary  processing  routines  into  files  (called 
B-files)  suitable  for  tape  storage. 

The  B-files  contain  only  essential  MESA  data  - temperature-corrected 
acceleration  counts,  derived  density  values  and  the  associated  GMT 
time  tags. 

It  should  be  noted  that  should  the  need  arise,  the  primary  routines  have 
the  capability  of  reading  this  B-ftle  as  input  to  perform  their  tasks. 
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Figure  15.  Printer  Report  Generated  by  the  GU-list  Section  of  MODPAP  (Cont. ) 
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This  capability  alleviates  the  need  to  repromote  TM  data  at  some  later 
date  should  a particular  orbit  of  data  require  reprocessing. 

2.1.8  FILM  — The  FILM  routine  provides  the  capability  of  displaying 
any  combination  of  the  following  35  mm  (or  16  mm)  film  plots: 

(a)  raw  MESA  data  displayed  as  a function  of  GMT, 

(b)  temperature-corrected  edited  MESA  data  as  function  of  GMT, 

(c)  drag  components  as  a function  of  GMT, 

(d)  atmospheric  density  as  a function  of  altitude, 

(e)  ratio  of  measured  density  to  model  density  as  a function  of 
altitude. 

Input  to  this  routine  is  the  GU  file  generated  by  RAWDATA  and  DRAG. 
Figures  6 and  9 (Section  1.3)  are  examples  of  raw  data  and  density  plots 
generated  by  this  routine. 
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2.2  MESA  Utilities 


In  addition  to  the  primary  routines  comprising  the  MESA  DRS,  there 
are  three  important  programs  which  are  being  utilized  daily  to  facili- 
tate the  processing  and  management  of  the  MESA  data  system.  A brief 
description  of  each  follows. 

2. 2. 1 DEMOTE  - Because  of  extremely  limited  on-line  storage, 
Geophysical  Unit  (GU)  files  must  be  demoted  (i.e. , copied)  to  magnetic 
tape.  For  obvious  reasons  it  was  desirable  to  group  processed  data 
orbits  by  satellite  and  subdivide  each  satellite  by  file  type  reduction 
status  on  separate  tapes.  For  this  reason  the  MESA  DEMOTE  facility 
was  developed. 

Simply  stated,  as  each  DRS  production  job  is  processed  the  final  step 
is  to  write  the  name  of  the  output  GU-file  into  a permanent  on-line  file 
called  DEMOTELIST.  The  DEMOTE  program  when  executed,  reads 
the  list  of  the  GU-files  to  be  demoted  from  DEMOTELIST.  Then, 
utilizing  information  acquired  from  the  Data  Management  Facility  tape 
directory  file  and  the  MESA  tape  selection  algorithm,  the  DEMOTE 
program  sets  up  input  files  containing  the  names  of  files  which  are  to 
be  demoted.  A listing  is  generated  describing  the  decisions  made  and 
a report  is  written  defining. the  tape  base  status. 

The  DEMOTE  processor  then  demotes  the  GU-files  to  the  tapes  chosen, 
and  the  original  DEMOTELIST  file  is  reinitialized  for  future  use.  This 
program  is  exercised  a few  times  a day,  and  it  provides  a reasonable 
method  of  maintaining  on-line  storage  space  for  MESA  processing. 


2. 2. 2 WORK  LIST  - The  WORKLIST  facility  is  a self-batching  proc- 
essor used  to  satisfy  most  MESA  DRS  requests  which  are  not  processed 
by  the  NASA  production  service.  It  provides  control  over  each  MESA 
DRS’  primary  processing  routine  and  GU-file  management  option,  as 


h 
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weJU  as  allows  the  user  to  directly  access  several  system  processors, 
WORKLIST  builds  a multitasking  environment  and  is  largely  self-main- 
taining, continuing  to  run  until  it  is  explicitly  stopped. 

Two  job  control  language  (JCL)  files  support  WORKLIST  by  batching 
invocations  of  each  other.  The  control  file  is  $:TASKS,  which  reads, 
interprets,  and  deletes  the  top  job  request  from  WORKLIST,  runs 
DRIVER9  to  prepare  NAMELISTS  for  the  MESA  primary  routines,  pro- 
motes any  GU-files  required,  batches  the  production  job,  $*TASK,  and 
finally  executes  any  processor  invoked  from  the  previous  job-set.  The 
production  job  maintains  the  WORKLIST  thread  by  rebatching  the  setup 
job,  $:TASKS  to  setup  the  next  request,  runs  each  routine  in  the  MESA 
DRS,  promoting  and  removing  telemetry  data  when  required,  runs  the 
various  processors  associated  with  GU-file  disposition,  and  finally  does 
clean-up  operations  to  remove  utility  files  associated  with  the  setup/run 
job-pair. 

To  ensure  that  output  files  from  separate  WORKLIST  jobs  have  different 
identifiers,  each  job  pair  utilizes  a unique,  2-digit  serial  number.  This 
number  is  updated  each  time  $:TASKS  is  run,  through  MESA's  ROTATE 
facility,  which  was  described  by  Section  2.1.1. 

The  use  of  WORKLIST  has  greatly  increased  the  processing  capability  of 
MESA  personnel,  while  at  the  same  time  requiring  minimal  manual 
maintenance. 

2. 2. 3 HEADHIST  — Due  to  the  large  data  reduction  requirements  of 
AE,  the  HEADHIST  facility  was  written  to  provide  an  automatic  logging 
system  for  all  data  as  it  is  being  processed  by  the  MESA  DRS.  HEADHIST 
is  a system  of  directory  files  containing  two  types  of  fixed  length  records, 
directory  records  a \d  header  records.  For  each  satellite  there  is  a 
separate  file  capable  of  storing  up  to  12, 000  orbits  of  directory  and 
header  information. 


I 


The  HEADHIST  system  is  accessed  through  calls  to  three  entry  points: 

(a)  HEADHIST,  which  stores  the  header  record  data  for  each 
orbit  processed, 

(b)  DIRECTORY,  which  updates  the  directory  records  as  each 
primary  routine  is  run, 

(c)  RESETDIR,  which  resets  directory  records  if  data  are  to 
be  removed. 

The  HEADHIST  files  are  read  by  other  programs  which  provide  a reason- 
able reporting  system  of  the  status  and/or  availability  of  reduced  MESA 
data.  Figure  17  is  an  example  of  one  such  report  generated  from  the 
HEADHIST  header  record  files.  An  example  of  other  programs  which 
access  the  HEADHIST  files  is  a program  which  generates  a listing  of 
calculated  bias  values  for  each  despun  orbit  processed  by  MESA.  These 
bias  data  will  be  later  used  to  study  the  performance  of  instrument  bias 
over  many  months. 


T I NG  HE  AL'ER  wATA  ’"*1  THL  A Il.f  rx^LtkLK-C  SPACECRAFT 


Figure  17.  Report  Generated  from  the  HEADHIST  Files 
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In  addition  to  the  MESA  DRS  and  MESA  utility  programs,  there  are 
other  MESA  programs  which  were  written  to  further  analyze  experi- 
mental data.  We  will  present  here  a brief  description  of  some  of  those 
used  for  data  analysis. 

2. 3. 1 ATTCQM  - In  order  to  evaluate  the  accuracy  of  the  spacecraft 
generated  attitude  data,  a routine  was  written  to  compare  M ESA-gen- 
erated attitude  data  with  the  spacecraft  generated  attitude.  MESA  atti- 
tude is  calculated  between  the  altitudes  of  180  - 200  km,  since  C-range 
MESA  data  saturates  at  lower  spacecraft  altitudes.  Algebraic  differ- 
ences are  calculated  for  this  region,  and  a printed  report  is  generated. 

This  routine  may  be  run  against  any  orbit  of  data  where  spacecraft  atti- 
tude values  are  suspect. 

2.3.2  NORMALS  - In  analyzing  density  data,  it  is  helpful  to  attempt 
to  remove  certain  effects  in  order  to  study  effects  due  to  other  causes. 
To  this  end  a computer  routine  was  written  which  would  normalize  cal- 
culated density  data  to  a specific  altitude  and/or  local  time  which  in 
effect  would  remove  density  effects  due  to  these  parameters.  Resultant 
normalized  density  may  then  be  more  easily  analysed  for  other  effects 
(such  as  magnetic  index  effects).  Normalized  density  is  calculated  by: 


measure  density  x model  density 

Normalized  density  = M the  normalized  altitude  and/o^local  time) 
J model  density 

Model  density  values  are  calculated  using  the  JACCHIA  71  Model 
atmosphere. 


Printer  listings  of  normalized  density  data  as  a function  of  time  and 
satellite  position  information  are  generated  from  this  program. 
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Figure  18  is  an  example  of  this  listing  showing  data  normalized  to 
an  altitude  of  200  km  and  a local  time  of  7. 5 hours. 

2. 3. 3 UAINTER  - The  Unified  Abstract  (UA)  file  data  base  was  organ- 

ized in  order  to  store  simultaneous  measurements  from  as  many  expe- 
riments as  possible.  For  this  reason  the  UA  file  was  setup  as  a function 
of  Greenwich  Mean  Time  (GMT)  at  specific  15  second  increments.  In 
order  to  compare  data  over  many  orbits,  however,  it  was  determined 
that  comparison  by  satellite  altitude  rather  than  by  GMT  should  be  done. 

Hence,  a computer  program  was  written  which  reads  data  from  both  the 
UA  point-by-point  files  and  the  spacecraft  orbit/attitude  (OA)  files. 

MESA  density  data  are  then  interpolated  at  specific  altitudes  for  each 
orbit  of  data  stored  in  the  UA  files.  Interpolated  MESA  density  values 
along  with  associated  orbit/attitude  position  data  were  punched  on  data 
cards  and  shipped  to  AFGL. 

These  data  cards  were  then  utilized  as  input  to  various  other  analysis 
programs  written  for  the  AFGL  CDC-6600  computer  system.  Examples 
of  this  are  given  in  Figures  19  and  20.  In  Figure  19  MESA  density 
values  at  140  - 200  km  are  displayed  along  with  magnetic  index  Kp  and 
satellite  perigee  latitudes.  In  Figure  20  MESA  density  data  for  140  km 
are  compared  to  Jacchia  71  model  values. 

In  addition  to  generating  MESA  density  data  cards,  UAINTER  was  modi- 
fied to  have  the  capability  of  extracting  from  the  UA  file  neutral  atmos- 
phere composition  data  calculated  by  the  OSS  experiment.  OSS  calcu- 
lated values  of  nitrogen,  total  oxygen,  and  argon  are  extracted  for  those 
orbits  where  MESA  density  values  are  stored.  This  will  provide  com- 
parison data  of  neutral  mass  density  from  both  experiments,  as  well  as 
determining  a way  to  provide  more  accurate  knowledge  of  the  lower 
thermosphere. 
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Figure  18.  Printer  Listing  Generated  by  NORMALS 
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Figure  19. 


AE-C  140  km  MESA  Density  Data,  Magnetic  Index  Kp,  and 
Perigee  Latitude 
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Figure  20.  AE-C  140  km  MESA  Density  Data  and  Jacchia  71  Model  Data 
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2. 3. 4 STATS  - An  alternate  method  of  separating  density  effects  due 
to  such  causes  as  Kp,  latitude,  longitude,  and  invariant  latitude  was  pro- 
grammed in  a routine  called  STATS.  This  program  utilized  the  data 
base  generated  from  UAINTER  (See  Section  2. 3. 3)  as  input,  and  it 
attempts  to  separate  density  points  into  particular  atmospheric  param- 
eter cells.  An  example  of  this  would  be  to  consider  all  density  points 
for  AE-C,  160  km,  between  the  following  values: 

60°  - 70°  latitude 

4-9  Kp 

-180°  - +180°  longitude 

0-12  hrs  GMT 

For  that  particular  set  of  parameters,  all  density  values  would  be 
analyzed,  their  average  value  would  be  calculated,  the  ratio  of  this 
average  value  to  an  average  value  of  the  Jacchia  71  model  (with  the 
same  parameters)  would  be  calculated  as  well  as  the  standard  devia- 
tion. 

Printer  reports  are  then  generated  showing  results.  Figure  21  illustrates 
this  type  of  output. 

STATS  is  a very  useful  tool  in  analyzing  these  data,  since  any  number 
of  different  parameter  cells  may  be  chosen. 

2.3. 5 NUT  AT  - Due  to  stringent  instrument  requirements  on  AE  and 
the  fact  that  the  AE  satellites  will  dip  into  very  low-perigee  orbits,  it 
is  essential  that  the  spacecraft  be  quite  stabilized,  that  is,  that  vehicle 
motions  such  as  nutation  be  kept  to  a minimum.  On  AE,  nutation  angle 
buildup  may  be  controlled  either  by  use  of  on-board  dampers  or  by 
judicious  use  of  momentum  wheel  speed,  hi  order  to  effect  these  con- 
trols it  is  necessary  to  monitor  their  results.  This  may  be  accomplished 
by  use  of  the  MESA  experiment. 
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Figure  21.  Observed/Model  Density  Statistics 
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We  developed  a printer  plotting  routine  to  display  nutation  amplitude 
and  frequency  as  a function  of  GMT,  by  displaying  the  MESA  XY-sensor 
output  data  taken  above  300  km.  From  this  preliminary  nutation  informa- 
tion may  be  determined  from  inspection  and  some  minor  calculations. 

However,  since  more  automated  techniques  were  desired,  we  developed 
a routine  called  NUTAT,  which  using  raw  MESA  sensor  data  as  input, 
calculated  nutation  angle  and  frequency  directly.  This  routine  utilized 
mathematical  modelling  and  optimization  techniques.  Results  agreed 
favorably  with  those  calculated  by  the  initial  printer  plot  routine. 


'1 


I 
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As  a result  of  this  contract,  we  have  processed  over  2000  orbits  of 
AE-C  MESA  telemetry  data  into  geophysical  parameters  and  stored 
these  results  in  files  (called  GU  files)  on  the  Sigma-9  computer.  In 
addition,  in  support  of  the  AE  Aeronomy  Team,  data  from  more  than 
1000  of  these  orbits  have  been  written  into  Unified  Abstract  (UA) 
files  which  are  accessible  for  use  by  all  other  experimenters  and 
theoretical  analysts. 


In  the  short  time  since  the  launches  of  AE-D  and  AE-E,  260  orbits 
of  AE-D  MESA  data  and  80  orbits  of  AE-E  data  have  been  processed 
and  stored  in  GU  files.  Of  these,  200  orbits  of  AE-D  data  and  65 
orbits  of  AE-E  data  have  been  written  into  UA  files  for  use  by  other 
team  members. 


We  teel  that  the  data  bases  generated  from  this  project  will  be  ex- 
tremely useful  in  any  further  studies  of  the  phenomena  of  the  upper 
atmosphere. 
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